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The opening of the pseudogap in underdoped 
cuprates breaks up the Fermi surface, which may 
lead to a breakup of the d-wave order parameter 
into two subband amplitudes and a low energy 
Leggett mode due to phase fluctuations between 
them. This causes a large increase in the tem¬ 
perature range of superconducting fluctuations 
with an overdamped Leggett mode. Almost res¬ 
onant scattering of inter-subband phonons to a 
state with a pair of Leggett modes causes anoma¬ 
lously strong damping. In the ordered state, the 
Leggett mode develops a finite energy, suppress¬ 
ing the anomalous phonon damping but leading 
to an anomaly in the phonon dispersion. 

The unexpected discovery of a giant anomaly in the 
dispersion of low energy phonons (GPA) in underdoped 
pseudogap cuprates has stimulated reconsideration of the 
role of phonons in cuprate high-Tc superconductors 
Recently many groups have proposed these anomalies are 
caused by other electronic instabilities e.g. charge density 
wave (CDW) order and also pair density wave (PDW) 
order, which compete with the uniform d-wave pairing 
state |6HI3]. A novel proposal has been put forward by 
Lee, who argues that Amperean pairing occurs in the 
pseudogap phase leading to an instability towards PDW 
and also CDW order mi- 

Two recent studies of YBCO samples covering a range 
of hole densities, found an onset hole density 0.18 

for the lattice anomalies, which coincides with the onset 
of the pseudogap usms]. Early photoemission (ARPES) 
experiments found that the onset of the pseudogap is 
characterized by a breakup of the Fermi surface into 4 
pieces centered on the nodal directions m- A rapid 
expansion of the temperature range of superconducting 
fluctuations above the transition temperature for long 
range superconductivity, Tc{p), is also observed |18| . 
This combination of the onset in hole doping, coinciding 
with Fermi surface breakup, and the onset in tempera¬ 
ture, coinciding with the onset of SC fluctuations, leads 
us to examine possible consequences of the special dis¬ 
connected nature of the Fermi surface in the pseudogap 
phase, on d-wave superconductivity. We And that su¬ 
perconducting fluctuations in an extended temperature 
above Tc can result as a special feature of d-wave su¬ 
perconductivity in the presence of the pseudogap. We 


shall show below that these fluctuations in turn can cou¬ 
ple to finite wavevector phonons leading to CPA. In our 
model the lattice fluctuations are dynamic as argued by 
LeTacon et al HEl] but random static CDW can still 
be induced by local perturbations, such as the random 
acceptors in nearly all underdoped cuprates |16| . The 
stoichiometric underdoped cuprate YBa 2 Cu 408 is an ex¬ 
ception. The NMR/NQR experiments by Suter et al [19] 
found dynamic charge fluctuations but no static lattice 
ordered modulation in agreement with earlier NMR ex¬ 
periments [201 nT\ . A recent detailed NMR study of 
YBCO found evidence for static lattice distortions, pos¬ 
sibly induced around lattice imperfections by CPA. They 
also did not report systematic splitting of the NMR 
lines which would be evidence for long range ordered 
CDW US]. 

In our microscopic model the key ingredient is the novel 
breakup of the full Fermi surface into 4 disconnected 
pieces that characterizes the pseudogap phase. In the su¬ 
perconducting phase, the 4 pieces combine to give 2 sets 
of Cooper pairs with nodes along (1,1) and (1,1). This 
breakup opens the possibility of superconducting phase 
fluctuations not just of the overall Josephson phase, but 
between the disconnected sets of Cooper pairs. The pos¬ 
sibility of phase fluctuations between separated pieces of 
the Fermi surface in multiband s-wave superconductors, 
i.e. the Leggett mode (LM), was studied by Leggett |22] 
many years ago. The superconductivity in cuprates has 
d-wave symmetry with nodes on the pockets and the 
accompanying sign changes reduce the net coupling be¬ 
tween the two sets. We propose this can lead to a low 
energy LM. 


Representation of Fermi surface 

One-loop Renormalization Croup calculations of the 
single band Hubbard model near 1 / 2 -fllling, show that 
at low energies the Coulomb repulsion develops substan¬ 
tial structure in k-space, peaking at momentum transfers 
(tt, tt) in both the particle-hole spin triplet and particle- 
particle singlet channels IM1I21]. The latter drives d-wave 
pairing with maximum amplitudes at antinodal. Yang 
et al |25l IM] argued the transition into the pseudogap 
phase is driven by increasing Umklapp scattering, which 
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Figure 1: Representation of the band structnre by 
Fermi arcs, a, The breakup of Fermi surface to 2 subbands 
a & 6 in (1,1) & (1,1) directions. Qgpa is a wavevector con¬ 
necting the two subbands, which is also the wavevector of 
phonon anomaly, b, Simplified model of Fermi arcs. Each 
Fermi arc is represented by a circular arc (shown red) with 
center Y, radius R, and terminates at the Umklapp surface. 
The Fermi velocity vf (blue arrow) is assumed to have con¬ 
stant magnitude on the whole arc. Y is uniquely determined 
by the choices of the wavevector between arc tips to be O.SItt, 
and the hole concentration p = 11.5%. These are typical val¬ 
ues in ARPES experiment. ± denotes the sign of symmetry 
factor 7 p. 


converts the antinodal gaps to insulating. They put for¬ 
ward an ansatz based on a modified BCS self energy with 
an energy gap pinned to the Umklapp surface leading 
to anisotropic Fermi pockets centered on the nodal di¬ 
rections. Later detailed ARPES experiments confirmed 
anisotropic nodal Fermi pockets m- The superconduct¬ 
ing complex pairing amplitude is confined to 2 discon¬ 
nected pairs of pockets centered on the (1,1) & (1,1) 
directions, illustrated in Fig. [T] We shall refer to these as 
subband a & 6 respectively. An examination of the pair 
scattering processes shows that there are intra-subband 
(tTjTt) processes. But the inter-subband processes con¬ 
necting subbands a & 6 involve cancellations between 
pairing and depairing scattering processes with similar 
wavevectors in a d-wave state. For this reason when we 
write the pairing scattering in terms of intra-subband 
and inter-subband processes for a d-wave state, we pro¬ 
pose that the intra-subband pairing is stronger than the 
inter-subband pairing, similar to the pairing model ex¬ 
amined by Leggett |22| . 

For simplicity the anisotropic nodal Fermi pockets will 
be represented as Fermi arcs (see Fig.[^ with a constant 
Fermi velocity along the arc. The parameters are chosen 
from a recent paper by Comin et al |28| . 


Leggett Mode and Fluctuations 


Assuming that both intra-subband and inter-subband 
couplings U and J are separable d-wave forms with sym¬ 
metry factor, 7 p, allows us to obtain the fluctuation pair 


propagator from the following Bethe-Salpeter equation 
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for Lq, in terms of which the full anisotropic fluctuation is 
written as Lpp>q = jpLqjpi with q = (q, fgo)- The elec¬ 
tronic bubble 4*^ = 

is defined on both sets of Fermi arcs for i = a,b. We de¬ 
note zero-temperature, finite-temperature, and retarded 
Green’s functions by Gp ^^, Gp^iu,, and Gp ^ respectively, 
similarly for other quantities. In the temperature region 
Tc < T < To (onset of superconducting fluctuations), 
the electron’s retarded Green’s function takes the form 
G?,. = [ w — Ep -I- , where ep is the Fermi arc dis¬ 

persion and F*^®) = aT -f hT^ is a temperature dependent 
quasiparticle damping [23]. The retarded propagator for 
the fluctuating LM is derived to be 




( 2 ) 


See the supplementary information (SI) for details. 
This form of the fluctuation propagator describes an 
overdamped bosonic mode. Here Ag is the density 
of states per spin at the Fermi energy for one pair 
of arcs, the constant c = W T ^ {a + bTc)\ 

with '0 (x) the digamma function. The damping of 
the LM = T~^ + Dq^, in which the inverse 

relaxation time ^ log -I- 2bT {T — To) + 

^^2 I J| / — J^) and the diffusion constant 

^ [h ++ /'P'[h ++ 

are functions of the temperature. The uniform damp¬ 
ing rate decreases as long range order at Tc is 
approached. 

We take zero temperature as a representative 
case for the ordered phase. In this case the 
LM is derived similarly to equation Q, but with 
a different form of the electronic bubble pa Ea 



^^+p,s^+ui^^-p,^-uiy ^PMi ~ T Cp) /Z Fp^^ — 
Ap/Z are the normal and anomalous Green’s functions, 
where Z = uj'^ — Ep+i5 and 5 —>• 0+. Ep = ^ 

is the quasiparticle energy and Ap = 7 pA is the d-wave 
gap function. It follows (see SI) 
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with a LM dispersion Wq = Wg + ^v^q'^. In this case, 
the LM is a coherent bosonic mode with infinite lifetime. 
The LM is gapped, and its frequency at zero momentum 
satisfies uJq = ^j-2|J|/ ([/^ —J^). Usually the ratio 

2A ^ 
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Figure 2: Feynman diagram of phonon self energy dne 
to the interaction with Leggett mode. B is the effective 
interaction vertex consisting of a particle-particle electronic 
bubble (G is the electron Green’s function). The intermediate 
state consists of two Leggett modes (L) and one phonon (D). 
Due to the separation of energy scales between Leggett mode 
and the electron, all quantities in green are neglected. 


Phonon Self Energy 


The k-space separation of the two bands does not move 
the LM away from q = 0, since this phase mode involves 
the transfer of zero momentum Cooper pairs between the 
subbands. Nonetheless it involves moving charges be¬ 
tween the subbands. Absorption and emission of phonons 
with the appropriate wavevectors also causes a charge 
transfer, but now as single quasiparticles, between the 
two subbands. Therefore it is not unexpected that a cou¬ 
pling between these processes should exist. In particular, 
we find the coupling is largest in the temperature region 
T ~ Tc, where the LM drops to zero energy and becomes 
overdamped. To this end we consider the process out¬ 
lined in Fig. where an incoming phonon is scattered to 
a nearby phonon wavevector with emission and absorp¬ 
tion of LM fluctuations. Such a process does not occur 
in standard superconductors but can exist here because 
of a soft overdamped LM for < T < Tg. Below we 
summarize calculations of the phonon self energy in two 
temperature regions. 

First we look at the phonon damping in the range of 
strong SC fluctuations starting at the onset temperature 
To of the SC fluctuations down to the superconduct¬ 
ing transition temperature T),. The expression for the 
phonon self energy 11, corresponding to the Feynman di¬ 
agram Fig.|^ follows 
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where aq is the electron-phonon coupling constant, 
Iq.k.iQ is the frequency summation over the intermedi¬ 
ate state consisting of two LMs and one phonon, Bq ^q 
is the effective interaction vertex between phonons and 
LMs, and DiQ = 2no/ — fig is the bare Green’s 

function for phonons with an assumed flat dispersion 
riq = rig. We have chosen the simplest form of the effec¬ 
tive interaction, B and ignore damping due to quasipar¬ 
ticle excitations, keeping only ReSq^, to concentrate 
on the novel phonon damping caused by the presence 
of a soft LM. Note, ReRq has a strong dependence 
on the phonon wavevector Q (see Fig. and peaks at 
a wavevector joining the ends of the arcs, because the 
symmetry factors and the available phase space for the 
transition at this wavevector are both large. We checked 
that ImRQ q <C Rei?o n parameters. 


ReRq also shows a peak at Q = 0 which will be dis¬ 
cussed later. 

The effective interaction vertex B involves an inte¬ 
gration over the whole Brillouin zone and all frequen¬ 
cies, while the LM L is only well-defined for small mo¬ 
menta and frequencies. This leads to a separation of 
spatial and temporal scales and enables us to ignore 
all small wavevectors and frequencies (marked green in 
Fig.§ in calculating B. Aided by the similarity with the 
Aslamazov-Larkin diagram |32fl36] . we perform analyti¬ 
cal calculation for Jq . For |fl — r2g| ^ fig, it follows (see 
SI for details) 
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In this region the on-shell values satisfy Re/jf^ = 0 and 
ImJ^^ < 0. The temperature dependence of the imagi¬ 
nary part of the retarded phonon self energy Imll^ i.e. 
the phonon damping, is plotted in Fig.|^ The self energy 
has a peak in momentum space at Q = (QgpAjO), near 
to the tip to tip wavevector between two sets of Fermi 
arcs. Because of the factor r in equation (|^, the tem¬ 
perature dependence shows anomalous behavior at the 
long range critical temperature Tc, in agreement with 
experiment |1]. 
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Figure 3: GPA in the fluctuation region, a, Temperature 
and momentum dependence of the effective interaction vertex 
B, with the dashed line marking Qgpa- b, The anomalous 
phonon damping. Parameters: Fermi velocity vf = SOOmeV, 
bare phonon frequency flo = lOmeV, ratio between Leggett 
mode frequency at q = 0 and the superconducting gap ^ = 
0.1, quasiparticle damping F^®^ = 0.5T + (0.3meV“^) and 
the long range order temperature Tc = 50K. We used an 
energy cutoff of TlOOmeV around the Fermi surface, which 
does not affect the qualitative feature of GPA. 


Below Tc, there is a finite restoring force for inter¬ 
subband phase fluctuations and the LM develops a finite 
energy at q = 0, which raises the energy of the inter¬ 
mediate state in Fig. As a consequence the approx¬ 
imate resonant condition between the incoming phonon 
and the intermediate state with a scattered phonon and 
2 LMs no longer holds, leading to a suppression of the 
phonon damping at low T. The GPA changes its form at 
T < Tc with strongly reduced damping. An anomaly in 
the phonon dispersion appears due to virtual coupling to 
an excited intermediate state. Treating the low tempera¬ 
ture behavior at T = 0, the factor from the intermediate 
state becomes 


Summary and Conclusions 


As we remarked earlier the transition into the pseudo¬ 
gap phase at the hole density pd ~ 0.18 displays three 
strong anomalies simultaneously — an antinodal insulat¬ 
ing energy gap leading to a breakup of the Fermi surface 
into 4 nodal pockets, a rapid expansion of the tempera¬ 
ture range of superconducting fluctuations, and the ap¬ 
pearance of a giant phonon anomaly in this temperature 
range. These phenomena are unique to the underdoped 
cuprates. Our aim here is to put forward a microscopic 
scenario which explains the interrelation between these 
phenomena. The low energy Coulomb interaction must 
be strong to drive the partial truncation of the Fermi 
surface, which is a precursor to a full Mott gap at zero 
doping |25l . A weakness of our microscopic scenario 

is the need to assume a form for this effective Coulomb 
interaction. Our choice is guided by the evolution found 
by Functional Renormalization Group in the overdoped 
density region |231 [21] . The persistence of d-wave sym¬ 
metry even as the maximally gapped antinodal regions 
transform from a superconducting to an insulating gap, 
leads to the conditions for a low energy LM to emerge. 
As we discuss above, this novel assumption enables us 
to consistently explain all three anomalies. In particu¬ 
lar we can explain the temperature evolution of the GPA 
from increasing damping as T ^ Tc from above, to a 
GPA with vanishing damping but a dispersion anomaly 
at T < Tc- Here we considered only zero magnetic field. 
The recent analyses of quantum oscillation experiments 
at high magnetic fields are consistently explained by a 
coherent orbit around all four arcs, is intriguing |3Z|. It 
raises the question of the evolution of the LM with in¬ 
creasing magnetic field for future study. 
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where Dq = 2D,o / (fP — LIq -b i6). We conclude that in 
this region, the on-shell values satisfy Re/jf^j < 0 and 
= 0, in agreement with the experiment for T < 

Tc 0. 

At long wavelengths, Q ~ 0, the long range nature of 
the Coulomb interaction suppresses the response in met¬ 
als at low frequencies to any perturbation coupling to the 
total electronic density. The electron-phonon interaction 
introduced in equation couples equally to both sub¬ 
bands and as a result the associated scattering processes 
are suppressed at Q ~ 0. 
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